Blockade and Counterflow Supercurrent in exciton-condensate Josephson junctions 
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We demonstrate that perfect conversion between charged supercurrents in superconductors and 
neutral supercurrents in electron-hole pair condensates is possible via a new Andreev-like scattering 
mechanism. As a result, when two superconducting circuits are coupled through a bilayer exciton 
condensate, the superflow in both layers is drastically modified. Depending on the phase biases the 
supercurrents can be completely blocked or exhibit perfect drag. 



Introduction — The terms superconductivity and su- 
perfluidity refer to dissipationless flow in charged and 
neutral systems, respectively. Superfluid exciton con- 
densates, in which macroscopic phase coherence is es- 
tablished among pairs composed of electrons and holes 
in different bands [U [2] , have been realized only recently. 
Signatures of exciton condensation have been reported in 
quantum Hall bilayers in which electrons and holes 
are located in two separate two-dimensional electron lay- 
ers [1], and in optically-excited exciton [5] and exciton- 
polariton [Hj cold gases. When the two layers of a bilayer 
exciton condensate (EC) are contacted separately, it can 
exhibit remarkable transport anomalies O [71 [8] associ- 
ated with its neutral supercurrents [3]. These properties 
provide an appealing platform for spectacular electrical 
effects in EC-superconductor hybrid systems in which 
the charged superconducting order parameter interfaces 
with the neutral EC order parameter. In this Letter we 
demonstrate that when two superconducting circuits are 
coupled through a bilayer EC, the superflow in both lay- 
ers is drastically altered. If the same phase bias is applied 
to both junctions, no Josephson current can flow through 
the system, a phenomenon we refer to as exciton block- 
ade. When a phase bias is applied to only one layer, on 
the other hand, it induces a superdrag counterflow super- 
current of the same magnitude in the unbiased layer. 

In order to explore the physics of conversion be- 
tween EC and Cooper-pair supercurrents, we consider 
the superconductor-EC-superconductor (S-EC-S) setup 
sketched in Fig. [T] Two closely-spaced layers, assumed 
to host an EC, are independently contacted to four su- 
perconducting electrodes. The electrodes in each layer 
are separated by a distance L much larger than the ex- 
citon coherence length, and an independent phase bias 
is applied to the top and bottom contacts. In the pres- 
ence of these biases, Josephson currents flow through the 
double layer. Because the EC is gapped, only dissipation- 
less counterflow can contribute to the Josephson current 
when L is long. The EC and the dissipationless nature 
of its counterflow supercurrent can therefore be revealed 
by a purely coherent equilibrium measurement when con- 
tacted by superconducting electrodes. 

Exciton blockade and superdrag — All important fea- 
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FIG. 1. (Color on line) Sketch of a superconductor-exciton- 
condensate-superconductor system. Panel a) A double-layer 
exciton condensate is contacted with four superconducting 
leads. Panel b) When the applied phase biases in the top and 
bottom layers are identical {ip in this cartoon) no current can 
flow when the length of the junction L is much larger than 
the exciton-condensate coherence length /ii;p/|r|. In this case 
the Josephson currents are in the "exciton blockade" regime. 
Panel c) When a phase bias Lp is applied to the top layer 
only, a supercurrent / flows. In the limit L 3> fiup/jFj a 
supercurrent — / is dragged in the bottom layer in a perfectly 
frictionless manner: in this case one has a "perfect" drag or 
"superdrag" of Josephson currents. 



tures of the physics we want to describe are captured by 
the simple one-dimensional (ID) model that we now dis- 
cuss in explicit detail. Later we provide arguments sup- 
porting the general validity of the conclusions we reach. 
The four superconducting leads in Fig. [l] are in an equi- 
librium configuration, characterized by the same Fermi 
level Ep. In each layer the electron filling is controlled 
by a gate voltage. In our calculations we assume that 
the two layers are oppositely gated, so that one layer 
(the top layer, say) is "hole-doped" and the other (the 
bottom layer) is "electron-doped" , as depicted in Fig. [2] 
Notice that states near the right Fermi point -\-kp are 
right-movers (left-movers) for the bottom (top) layer, 
whereas states near the left Fermi point —kp are left- 
movers (right-movers) for the bottom (top) layer. Let 
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^aa be the field operator describing an electron in the 
a = top (T)/bottom (B) = ± layer with spin a =t, i- We 
have applied the Bogoliubov-de Gennes approach [TOl [11] 
to model our S-EC-S hybrid system: the Hamiltonian T-L 
of the system is ?^ = dx 'i/^ {x)'H{x)'i> (x) , where en- 
ergy is measured with respect to the equilibrium Fermi 
level Ef- Here ^ = and is a 

4x4 matrix 
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It contains single-particle band-kinetic-energy terms for 
each layer, intra-layer terms containing the supercon- 
ducting order parameter [TU] Aq, c>c (^'a^^a^) and inter- 
layer terms containing the EC order parameter [3] F cx 
(^Bcr^Tcr)- Both order parameters vary spatially along 
the current-flow {x) direction. Since details of order pa- 
rameter behavior near the interfaces are irrelevant for 
our purposes, we can assume that are non-vanishing 
only in the left (L) and right (R) electrodes. For sim- 
plicity we also assume the same amplitude |A| in all the 
electrodes, whereas the phases (/?q,l/r are allowed to dif- 
fer. More explicitly, we take A^ix) = | A|e''^° Le(-a;) -I- 
I Aje*"^" " 9(a; — L), where the origin a; = is chosen at 
the left layer-electrode interfaces and 8 is the Heaviside 
step function. In contrast, the EC order parameter is 
taken non- vanishing in the double-layer region, i.e. for 
Q < X < L. Although its amplitude |F| can be taken 
as constant, it is essential to allow for phase variation in 
order to account for condensate counterflow currents. In 
ID current conservation implies linear phase variation so 
that F has the form 



(2) 



When phase biases are applied to the four electrodes, 
supercurrents flow in both layers. The EC weak-link sup- 
ports two contributions to the Josephson current. The 
quasiparticle channel contribution, in which Cooper pairs 
propagate by the virtual excitation of quasiparticles in 
the double layer, is present in ordinary weak links. In 
the present case, however, it is exponentially suppressed 
when L ^ /ivp/lrj {vp being the Fermi velocity) because 
of the gap in the quasiparticle excitation spectrum of the 
EC. Much more interesting is the new contribution to the 
current which derives from the conversion of supercurrent 
into superfluid excitonic current. It can be visualized as 
a correlated Andreev reflection [12\ in which an electron 
and hole (in different layers) enter the EC and propa- 
gate without dissipation to the other end of the double 
layer. There a similar process occurs to convert the exci- 
ton current back into a Cooper-pair current. This process 
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FIG. 2. (Color on line) Schematic electronic band structure 
of a gated bilayer in the absence of electron-hole coupling. 



survives also in the long-junction limit and it leads to a 
number of spectacular effects, as we shall show below. 

In the long junction limit, |A|, |F| 3> hvp/L, the crit- 
ical current does not depend on the magnitude of ei- 
ther order parameter. The mathematical description 
of the long-junction limit is simplified if we also as- 
sume that |A| ^ |F|, with no physically relevant con- 
sequences for the main results. Indeed, for energies 
much smaller than |A| the layer spectra can be lin- 
earized [13, , and the electron field operators can be writ- 
ten as "i/aaix) = e''=''^^'a<T+(a;) + e' "^a.y-ix), where 
"^i^^j-i^x) are slowly varying fields related to the Fermi 
points ifcp. Furthermore, the presence of the supercon- 
ductors can be accounted for by boundary conditions at 
the contacts, such as {0} = ±ai e*"^"-^ ^l^(-t-)- (0) 

at the left interfaces. Similarly for the right interfaces lo- 
cated at X — L. In this way the problem is reduced 
to the evaluation of the average of the current operator 
over the equilibrium state of a system in which the fields 
satisfy these boundary conditions and the exciton order 
parameter ^ exhibits a space-dependent phase winding. 
It follows that the EC winding wavevector q must satisfy 
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with J an integer. Here Lpa = ^Pcch — Va,R is the phase 
bias in layer a. Minimization of the total energy fixes 
J to be the closest integer to {ipT ~ ¥'B)/(27r), and the 
offset phase to be 70 = {'Pt,l — ¥'bx)/2. We find that, at 
zero temperature, the supercurrents in top and bottom 
layers fiow in opposite directions. Explicitly, they exhibit 
a sawtooth form. 
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where ipT — ipB is defined modulo 27r. The magnitude of 
the currents depends on the difference ip^ — (/Sb between 
the phase biases ip^ and ipB in the two layers. 

Eq. ^ is the main result of this Letter and has several 
interesting physical implications: i) ip^ = (pB (parallel 
fiow). When the same phase biases are applied to the 
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two junctions no supercurrents can flow through the EC. 
In this case the Josephson currents experience an "exci- 
ton blockade"; ii) (pT = ^'Pb (counterflow). In this case 
the Josephson current flowing through the EC is maxi- 
mal, with a critical value equal to the critical current of 
a ballistic one-channel superconductor - normal metal - 
superconductor (S-N-S) junction; iii) ip^ — p and ipB = 
(superdrag). When current flows in one layer due to a 
phase bias in that layer, a current equal in magnitude 
but opposite in direction flows in the other layer. This 
is a consequence of perfect conversion of exciton current 
into supercurrent. Eq. Q can then be seen as a perfect 
drag effect for the supercurrent. 

The existence of a dissipationless (counterflow) chan- 
nel also has a spectacular impact on the temperature 
dependence of the critical current. Indeed we notice 
that the ground-state current Q has the same length 
dependence as that in a S-N-S junction [M]. At finite 
temperature it is possible to show that, in the regime 
hvp/L < k^T |r|, the critical current in the S-EC-S 
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where /? — l/{kBT) and Lth — fiv-p / (k-QT) is the thermal 
length. Note that the second term in square brackets in 
Eq. ([5]) is oc exp(— /3|r|): thus, as long as thermal fluctu- 
ations are dominated by the excitonic gap, the ground- 
state current is essentially unaffected by thermal fluctu- 
ations. This is another important result of this Letter. 
Notice that this occurs even when the thermal length 
ith is smaller than the length L of the junction. This 
is in striking contrast with the case of a S-N-S junction 
(or with the case of two decoupled layers), where the 
critical current /c is exponentially suppressed [T^, i.e. 
Ic oc exp(— L/Lth) for Lth L, due to thermal deco- 
herence affecting a single Andreev-reflection process. In 
the presence of the EC, Andreev processes coherently 
occurring in the two layers transform Cooper pairs into 
electron-hole pairs of the EC, which are protected from 
thermal decoherence by the excitonic gap. Thus in the 
temperature window Hvp/L ^ fceT" ^ |r| the EC coun- 
terflow channel is responsible for an exponential enhance- 
ment of the critical current with respect to the ordinary 
S-N-S case mentioned just above. 

Discussion — The only crucial assumption we made 
in the previous derivation is that the length L of the 
junction is much larger than the EC coherence length 
/iUF/|r|. For these reasons, the physical results obtained 
are not restricted to the specific ID model discussed 
above. The dependence of the Josephson current on the 
difference — 'Pb can be deduced from quite general 
arguments. Current conservation indeed implies that the 
supercurrent can be evaluated in the bulk of the layers, 
where the supercurrent is purely carried by the EC, pro- 



vided that the junction is long enough (L » ?iVF/|r|). 
We also emphasize that, due to the charge neutrality 
of the EC order parameter F, the supercurrents in the 
two layers are equal in magnitude and opposite in sign. 
In particular, the current is proportional to g, i.e. the 
phase winding in Eq. ([2|. The evaluation of the cur- 
rent thus reduces to the determination of the depen- 
dence of q on the superconducting phase biases (^t and 
tf-Q. Let us assume that the system is translationally- 
invariant in the transverse direction y, so that the or- 
der parameters /S.a{r) and F(r) depend on the longitu- 
dinal x-direction only. By applying the transformation 

*a<T(r) = exp {i[v?a,L + Wa,K " (Pa.l?)X / L] /2}'i! aa{r) , 

the superconducting phase biases (pa,i./K can be gauged 
away. The price to pay is twofold. Firstly, an effective 
vector potential = {^Qipa/L) x appears in each intra- 
layer Hamiltonian (here $0 = h/2e is the quantum of 
flux, associated with elementary charge 2e of the super- 
conducting order parameter A^). Secondly, the EC order 
parameter describing the inter-layer coupling transforms 
into 
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By observing that the energy scales characterizing the 
intra-layer and the inter-layer terms are hvp/L and |r|, 
respectively, it is straightforward to realize that for a 
long junction (L ^ /iup/|r|) the inter-layer terms play 
the major role in determining the equilibrium configura- 
tion. Energy minimization implies that the argument of 
the exponent proportional to x in Q vanishes, since in 
the new gauge the system is effectively phase unbiased. 
This fixes the winding q to be the one defined in Eq. (|3| , 
and proves that the currents just depend on ipj^ — (p-Q. 
The current-phase relationship is always of the sawtooth 
form. We stress that the above argument does not depend 
on the details of the experimental setup. In particular, 
it applies independently of the specific (parabolic or lin- 
ear) energy-momentum dispersion relation of the intra- 
layer kinetic Hamiltonian. Furthermore, it also holds if 
the contacts with the superconducting electrodes are not 
ideal and when their transparencies are different in the 
top and bottom layers. In the specific case of 2D layers 
with width W and highly-transparent contacts, the cur- 
rent is still given by Eq. Q provided that it is multiplied 
by the number k^W/A of transverse channels [TIHI6) . 

The unique properties of the conversion of EC cur- 
rents into charged supercurrents can be exploited for a 
number of possible applications. As an example we dis- 
cuss a configuration realized by closing the two supercon- 
ducting electrodes contacted to the (say) top layer into 
a ring-shaped rf-SQUID geometry, so that the phase dif- 
ference iy9T is directly related to the magnetic fiux <I>t by 
(/?T = 27r$T/'i?o + 27rn. In response to a magnetic field, 
an induced Josephson current /x fiows in the top layer 
and, according to Eq. Q, an opposite current /b = —It 
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flows in the bottom layer. Whenever the magnetic field 
changes the flux by a fluxon, the currents in both lay- 
ers are reversed. Current sign switches detected in the 
bottom layer count the fluxons present in the top layer 
ring. If the magnetic flux is generated by a monotonic 
analog input signal, the system effectively converts it into 
a sum of current switch pulses, i.e. to a digital signal. 
The system is therefore an analog-to-digital converter. A 
generalization to non-monotonic input signals can easily 
be achieved by using two double junctions. 

One important obstacle which presently stands in the 
way of observing these effects is the fact that equilib- 
rium exciton condensation has so far been observed only 
in Quantum Hall (QH) bilayers at total filling factor 
i^T = 1- QH systems necessarily have current-carrying 
gapless channels localized at their edges. In a QH bar 
geometry the edge channels will alter the physics we dis- 
cuss. Spontaneous coherence between conduction and 
valence band electrons in different semiconductor quan- 
tum well layers is however also expected [ITJ to occur 
at zero magnetic field when inter-layer interactions are 
strong. There are hints that the conditions necessary for 
coherence have been realized in some recent |18j semi- 
conductor bilayer experiments. Graphene bilayer sys- 
tems |19j are just starting to be examined for coher- 
ence effects and have a number of potentially impor- 
tant advantages, as pointed out recently by several re- 
searchers Because they are gapless and atom- 
ically 2D, the field-driven carrier densities that can be 
achieved are much larger than in the semiconductor case. 
Weaker dielectric screening and linearly dispersive con- 
duction and valence bands help to increase both interac- 
tion and disorder energy scales. Finally, graphene bands 
are nearly perfectly particle-hole symmetric, guarantee- 
ing the nearly perfect nesting between conduction and 
valence band Fermi surfaces which favors the coherent 
state. Progress in the realization of electrically 
isolated double-layer graphene sheets, either two-layers 
separated by a dielectric or rotated layers, is on-going. 
In the single-layer graphene case, it has already been 
demonstrated that it is possible to fabricate trans- 
parent interfaces between graphene and superconducting 
electrodes. In view of this progress our predictions are 
likely to be within experimental reach soon. 
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